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O ■ ABSTRACT 

o : 

^ Using sensitive centimeter-wave receivers mounted on the Owens Valley Radio Observatory and 

04 ' Berkeley-Illinois-Maryland- Association millimeter arrays, we have obtained interferometric measure- 

O , ments of the Sunyaev-Zel'dovich (SZ) effect toward massive galaxy clusters. We use the SZ data to 

determine the pressure distribution of the cluster gas and, in combination with published X-ray tem- 
peratures, to infer the gas mass and total gravitational mass of 18 clusters. The gas mass fraction, /g. 
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\l , is calculated for each cluster, and is extrapolated to the fiducial radius rsoo using the results of numer- 

ical simulations. The mean fg within r^QQ is 0.08llQQ5i/ij~oo (statistical uncertainty at 68% confidence 
^ ■ level, assuming i7Af=0.3, VLj^=Q.7). We discuss possible sources of systematic errors in the mean fg 

, measurement. 

' We derive an upper limit for from this sample under the assumption that the mass composition 

■ of clusters within rsoo reflects the universal mass composition: VLMh < fls/fg- The gas mass fractions 

depend on cosmology through the angular diameter distance and the v^qq correction factors. For a flat 
universe (r2A= 1 - and /i = 0.7, we flnd the measured gas mass fractions are consistent with less 
than 0.40, at 68% confldence. Including estimates of the baryons contained in galaxies and the baryons 
which failed to become bound during the cluster formation process, we find VIm ~0.25. 
Subject headings: cosmic background radiation — cosmology: observations, galaxies: 
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2 ' clusters-techniques: interferometric 

to 

1. INTRODUCTION where f^s is the ratio of baryon mass density in the uni- 

r^^ , fi - u -4- -fu- ii verse to the critical mass density. The cluster gas mass 

Clusters or galaxies, by virtue oi being the largest known ^ . . 

■ ■ T J T- . ■ , , 1 ri 1 fraction measurement can be used within the Big Bang 

virialized objects, are important probes ot large scale ° ° 

^ , , , , jiii 1-1 ji Nucleosynthesis (BBN) paradigm to constrain a/, 

structure and can be used to test cosmological models. \ ; f b 

Rich clusters are extremely massive, ~ 10^^ M©, as indi- <^b/ fg- (2) 

cated by the presence of strongly gravitationally lensed „i i r r-. ■ ^ • i i n-n-KT i i i- i 

i he value oi \Ib is constrained by BBN calculations and 

background galaxies, the large velocity dispersion (> 1000 ^, , r^i i i r_Li i- i 

^ the measurements ot the abundances ot the light elements 

km s ) in the member galaxies, and the high measured ,,,, , , ^nr.r, ■ ^ ^r,r^r\ n 

' \ (Wagoner et al. 1967; Gopi et al. 1995) as well as measure- 

temperature (> 5 keV) of the ionized intracluster gas. The X r 1 • i • r ii 

^ ' ments ot the spatial anisotropics ot the cosmic microwave 

mass composition on cluster mass scales is expected to re- background (CMB) (White et al. 1994; Hu et al. 1997). 

fleet the universal mass composition (White et al. 1993; mi i • i • i ; • i ■ ;i 

ihe luminous baryons m clusters are mamly m the 
Evrard 1997) . Under the fair sample assumption, then, . ^ , ^ ,. i-,r^T.,\ mi • i ^ 

' . . . . . gaseous intracluster medium (IGM). ihe gas mass IS about 

the cluster gas mass fraction, which is a lower limit to the , ^ -.11 .1 ,1 • ;• n 1 

an order 01 magnitude larger than the mass m optically ob- 

cluster's baryon fraction, fs, should reflect the universal j 1 ^ 1 • ti7i ^ 1 ^1AAo^ c 

f w I served cluster galaxies, 6.17., White et al. (1993); lorman & 

baryon fraction. Jones (1982). Hence, the gas mass is not only a lower limit 

to the cluster's baryonic mass, it is a reasonable estimate 

fgl^fs— jY^' of it. Although observations suggest that galaxy groups 

^Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02139 
^Department of Astronomy &; Astrophysics, 5640 S. Ellis Ave., University of Chicago, Chicago, IL 60637 
3Dept. of Space Science, SD50, NASA Marshall Space Flight Center, Huntsville, AL 35812 
^Department of Physics, University of Alabama, Huntsville, AL 35899 
^Department of Physics, University of California, Berkeley, CA 94720 

1 



and low mass clusters may have lost gas due to preheat- 
ing or post-collapse energy input (David et al. 1995; Mohr 
et al. 1999; Ponman ot al. 1996), the gas mass fraction in 
massive clusters (Tg > 5 kcV) appears to be constant. The 
gas mass fraction in massive clusters then provides a lower 
limit to the cluster baryon fraction, fg < fs- 

The ICM is hot, with electron temperatures, T,,, from 
~5 to 15 keV; rarefied, with peak electron number densi- 
ties of Ue — 10~^ Cm~^; and cools slowly {tcool > t Hubble), 

mainly via thermal Bremsstrahlung in the X-ray band. 

The ICM also produces a spectral distortion of the CMB 
known as the Sunyaev-Zel'dovich effect. The ICM mass 
fraction may be calculated from either of these observables. 
In addition to providing measurements of this important 
parameter with independent techniques, the two methods 
are fundamentally different in that the SZ effect is directly 
proportional to the integrated density of the gas while the 
X-ray brightness is proportional to the integrated square 
of the density. 

The X-ray surface brightness is proportional to the emis- 
sion measure, oc J nlA{Te)dl, where the integration is 
along the line of sight. Under simplifying assumptions, 
the gas mass can be calculated from an X-ray image de- 
projection. Since the sound crossing time of the cluster 
gas is typically much less than the dynamical time, one 
may reasonably assume that, in the absence of a recent 
merger, the cluster gas is relaxed in the cluster's poten- 
tial. Hydrodynamic simulations support this notion, e.g., 
Evrard et al. (1996). Under the assumption that the gas is 
in hydrostatic equilibrium (HSE), supported only by ther- 
mal pressure, the total binding mass follows from the gas 
density and temperature distribution, the latter of which 
may be determined with X-ray spectra. Gas mass frac- 
tions have been measured with this technique out to clus- 
ter radii of 1 Mpc or more (White & Fabian 1995; David 
et al. 1995; Neumann & Bohringer 1997; Squires et al. 
1997; Mohr et al. 1999). In an X-ray flux-limited sam- 
ple of 45 clusters, Mohr et al. (1999) measure the mean 
cluster gas mass fraction within approximately the virial 
radius to be (0.0749 ±0.0005)/i-3/2. Here, and throughout 
the paper, we use Ho = lOO/i km s^^ Mpc^^. 

In this work, we calculate cluster gas mass fractions us- 
ing spatially resolved measurements of the SZ effect to de- 
termine the gas density profile. The SZ effect in clusters is 
a spectral distortion of the CMB radiation due to inverse- 
Compton scattering of the relatively cool CMB photons 
off hot ICM electrons (Sunyaev & Zel'dovich 1972). At 
frequencies less than ^ 218 GHz, the intensity of the 
CMB radiation is diminished as compared to the unscat- 
tered CMB, and the SZ eSect is manifested as a brightness 
temperature decrement towards the cluster. This decre- 
ment, ATsz/TcMB, has a magnitude proportional to the 



Compton y-parameter, i.e., the total number of scatterers, 
weighted by their associated temperature, 

I nS)Te{l)dl, (3) 

meC^ J 

where k is Boltzmann's constant, ctt is the Thomson scat- 
tering cross section, rUe is the electron mass. We extract 
the cluster's gas mass fraction from a deprojection of the 
SZ effect data in a method analogous to the described X- 
ray HSE analysis. 

In addition to providing an additional measurement of 
fg, we note several points of difference between the X- 
ray and SZ analyses. Should significant large-scale or 
spatially- varying clumping of the ICM be present, the SZ 
image deprojection may look quite different from the X- 
ray deprojection. Clumping at scales below the resolution 
of the X-ray and SZ images could also result in a differ- 
ence of < rig >-^/^ / < rie > in the inferred gas mass. 
Also, the emission from the cores of relaxed clusters may 
be dominated by cooling flows, which complicate the in- 
terpretation of the X-ray data and may bias the result 
strongly if not taken into account (Allen & Fabian 1998; 
Mohr et al. 1999). The direct relationship between the 
SZ effect and the gas density also permits a surface gas 
mass fraction to be measured without image deprojection 
by comparing the projected or "surface" gas mass from 
the SZ observation to a measurement of the surface total 
cluster mass, for instance from gravitational lensing mod- 
els (cf Grego et al. (2000)). Because lensing observations 
are not available for all the clusters in our sample, and 
because we are interested in the gas mass fraction within 
clearly defined cluster radii, for this work we calculate fg 
with the dcprojcction/HSE method only. 

In this paper, we present cluster gas mass fractions 
based on SZ measurements made in the years 1994-1998, 
and a discussion of the implications of these measurements 
for cluster physics and for cosmology. In previous papers 
(Carlstrom et al. 1996; Grego et al. 2000), we have de- 
scribed the instrument constructed expressly to make such 
measurements, and the reduction and calibration methods 
for the SZ measurements. We give a brief review of this 
and discuss the cluster selection and observations in Sec- 
tion 2. In Section 3, we describe the procedure for fitting 
the SZ data to models for the cluster gas and extract- 
ing cluster gas masses and gas mass fractions, including 
a discussion of the possible systematic uncertainties. In 
Section 4 we present the derived gas masses and gas mass 
fractions for the entire cluster sample, compare these re- 
sults to other gas mass fraction work, and discuss the limits 
this work places on Q.m and plans for future work. 

2. INSTRUMENT AND OBSERVATIONS 

2.1. An Overview of the OVRO and BIMA 30 GHz SZ 
Observations 
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We wish to take advantage of the characteristicaUy low 
noise of interferometer systems while retaining sensitiv- 
ity on the large angular scales subtended by clusters. To 
do this, we integrated centimeter-wave receivers built ex- 
pressly for this purpose into the millimeter-wavelength in- 
terferometer systems at the Owens Valley Radio Obser- 
vatory (OVRO) Millimeter Array and at the Bcrkclcy- 
Illinois-Maryland Association (BIMA) Millimeter Array. 
The angular scale sampled by an interferometer element is 
6 ^ X/B, where B is the projected baseline, or telescope 
spacing as seen by the source. At our 1 centimeter op- 
erating wavelength, the compact telescope configurations 
effectively sample the angular scales of clusters while the 
fluxes of any contaminating pointlike sources in the field 
are simultaneously monitored with the longer baselines, so 
their time variability is not a source of uncertainty. Oper- 
ating the millimeter systems at ~ 10 times lower frequency 
than the design frequency also provides for very good op- 
tical performance. Both arrays allow for the elements to 
be placed in a wide variety of configurations. 

The receivers, which operate from 26 to 36 GHz, are 
based on cryogenically-coolcd 4-stage HEMT amplifiers 
and achieve receiver temperatures of 12-20 K at 28.5 
GHz. Results from this system are reported in Carlstrom 
et al. (1996, 2000); Grego et al. (2000); Holzapfcl et al. 
(2000a,b); Patel et al. (2000); Reese et al. (2000). 

2.2. The Interferometric Arrays 

We observed with this system at the OVRO Millime- 
ter Array in the summers of 1994-1996 and 1998. At 
OVRO, the weather was adequate for observing about 80% 
of the time. The aperture efficiency at 28.5 GHz, ~ 0.75, 
was measured with holographic techniques. The contribu- 
tion of the antenna to the system temperature, including 
spillover, is ~ 12-15 K, as measured from sky dips. The 
array of six 10.4 meter telescopes (four telescopes in 1994) 
is two-dimensional, with baselines ranging from 14 to 240 
meters. A general description of the OVRO millimeter ar- 
ray is provided in (Padin et al. 1991). The continuum mea- 
surements are made with the dual-channel analog correla- 
tor, each channel having an input bandwidth of 1 GHz. In 
1994, the SZ observations were made using a single chan- 
nel centered at 28.7 GHz. After 1994, the observations 
were made in single-sideband mode using two 1 GHz chan- 
nels, centered at 28.5 and 30 GHz. At OVRO's latitude, 
sufficient u-v coverage can be obtained for sources with 
declinations greater than —10° when two or three differ- 
ent telescope configurations are used. The primary beams 
for each telescope are measured holographically and are 
quite similar. The full width at half maximum (FWHM) 
of the beams differ maximally by five percent, and can be 
approximated as Gaussian with a FWHM of 252". 



We used the same receiver system at the BIMA Mil- 
limeter Array in the summers of 1996, 1997, and 1998. In 
1996, we used the six receivers on six of the BIMA tele- 
scopes; three additional receivers were constructed to use 
a total of nine of the ten BIMA telescopes in the sum- 
mers of 1997 and 1998. At BIMA, the contribution to the 
system temperature from the antenna is minimal, ^ 6 K. 
The aperture efficiency at 30 GHz with our receivers is ~ 
0.70. The BIMA array is two-dimensional, with baselines 
ranging from as short as 7.5 meters and as long as 1 kilo- 
meter. A general description of the BIMA interferometer 
is given in (Welch et al. 1996). We operate the hybrid dig- 
ital correlator in wideband mode (mode 8 in the notation 
of Welch et al. (1996)) covering 800 MHz with 2-bit sam- 
pling. Adequate u-v coverage can be obtained for sources 
with declinations greater than about —10° with one or two 
telescope configurations. The primary beams for each tele- 
scope were measured holographically and are very similar. 
The FWHM of the beams differ maximally by ~3% and 
can be approximated by Gaussian with a 396" FWHM. 

2.3. Data Reduction and Calibration 

2.3.1. OVRO Reduction 

Our observing strategy maximized usable observing 
time on the cluster while also providing reliable instru- 
ment calibration. During times the cluster was observ- 
able with minimal shadowing, we interleaved twenty-four 
minute observations of the cluster with observations of a 
nearby bright radio source (the gain calibrator) to moni- 
tor the stability of the interferometer's phase and ampli- 
tude response. The cluster and gain calibrator observa- 
tions were taken in several short segments (four and one 
minute integrations, respectively) to minimize the effect 
of short-term instabilities on observing efficiency. Either 
a planet or a time-stable bright radio source was observed 
to provide the absolute fiux scale for the measurement. 
This flux scale is based on Mars; if Mars was at least 15° 
above the horizon during the cluster observation, it was 
observed. 

The data were edited according to several criteria. Data 
taken with a telescope which was blocked by another tele- 
scope (shadowed) are removed from the data set. We use 
a conservative shadowing limit; data are discarded if the 
projected baseline is less than 1.05 times the telescope di- 
ameter. Also removed are data taken during poor weather 
as evidenced by poor phase; stablity and data affected by 
anomalous jumps in the instrument's phase. Any cluster 
data not bracketed by calibrator observations are also re- 
moved. 

Data calibration proceeds in two steps, gain calibration 
and absolute flux calibration. A time series of the gain 
calibrator's amplitude and phase in each baseline is ex- 



3 



amined with the MMA data reduction package (Scoville 
et al. 1993). The instrument response during the cluster 
observations is interpolated from a fit to this time series. 
The amplitude response generally varied less than a per- 
cent over an observation of many hours. The average gains 
for each baseline were quite stable from day to day. In the 
1994 and 1995 observing seasons, however, the receivers 
responded to linearly polarized light. Since some of the 
gain calibrators are linearly polarized at the 5-10% level, 
the measured amplitude of such calibrators changes with 
parallactic angle as well as instrument response. Only a 
few of the cluster observations are significantly affected, 
since two different gain calibrators were often used for a 
single cluster and in no case were both noticeably polar- 
ized. For the affected clusters, the average calibrator flux 
is used and the amplitude gain is assumed to be constant. 
The instrumental phase response typically only drifts a few 
tens of degrees over the course of a cluster observation. 

The absolute flux scale is determined relative to Mars. 
Mars' brightness temperature is predicted using a radia- 
tive transfer model for the whole disk brightness temper- 
ature (Rudy 1987) for each day of observation. The in- 
trinsic uncertainty of this model is expected to be ~ 2.5% 
and the uncertainty from input parameter uncertainties is 
about 3%, and so we estimate the accuracy of our absolute 
flux scale to be 4% at 90% confidence. We calculate the 
brightness temperature at the center frequency of our ob- 
served band; the brightness temperature varies less than a 
percent over our bandpass. The solid angle Mars subtends 
at each observation is determined from the equatorial and 
polar diameters of Mars reported in the Astronomical Al- 
manac. 

Goldin et al. (1997) compared the Rudy model to a ther- 
mal model for Mars, and find even with substantial extrap- 
olations in wavelength, the two models predict brightness 
temperatures for Mars consistent with each other. We also 
compare the Mars brightness temperature predict(xl by the 
Rudy model to those derived by Mason et al. (1999) based 
on absolute flux measurements of Gas A. They find Mars' 
brightness temparature at 32 GHz to be 196.01y'| ^o^" 
the May 1998 epoch. In our observing scheme, we deter- 
mine the brightness temperature separately for each day. 
The brightness temperature predicted by the Rudy model 
at 32 GHz varies from 194 K to 203 K during the month 
of May 1998. The brightness temperature for Mars varies 
less than 0.2% between 28 and 30 GHz. These compar- 
isons suggest that our primary calibration is accurate and 
is consistent with the primary calibration used by other 
groups. 

The fluxes of a set of primary calibrators were deter- 
mined using the predicted Mars flux. Since the amplitude 
gain of the instrument is stable with respect to time and 



telescope elevation, the observations of these calibrators 
and Mars need not be contiguous. In the case our primary 
calibrator is never observable at the same time as Mars, 
we bootstrap the flux from another primary calibrator. 
Over each of the month-long observation seasons, no time 
variation of the gain calibrator fluxes was evident. 

2.3.2. BIMA Reduction 

At BIMA, we use an observing scheme similar to that 
used at OVRO, interleaving observations of the gain cali- 
brator and cluster. 

The reduction proceeds much like the OVRO reduction 
with additional editing and passband calibration. Spectral 
channels with low signal-to-noise ratio or with spurious 
interference are discarded. Also edited out are shadowed 
data, data taken with obviously incorrect or irregular sys- 
tem temperatures, data taken when the telescope tracked 
incorrectly, and data contaminated by local interference. 
(These errors are flagged online at OVRO.) The spectral 
response of the instrument is determined from a passband 
calibrator, and then the spectral channels are vector aver- 
aged into one wideband channel. The gain calibration is 
then performed. 

Absolute flux calibration at BIMA evolved between the 
1997 and 1998 seasons. For the 1997 data, each of the Mars 
observations were reduced in the method described above, 
and the resultant amplitude and phase are SELFCALed. 
The amplitude response for each of the 9 telescopes is de- 
termined using the flux of Mars from the Rudy model. The 
gains were very stable over the two months of observing 
time, with an r.m.s. antenna gain variation of 1.2% for 
all telescopes all summer. With the knowledge that the 
amplitude response is steady, in 1998 the gains were de- 
rived in the first week of the BIMA observations, and these 
gains are applied online. Mars was subsequently observed 
to monitor any gain variation. 

2.4. Cluster Selection and Observations 

We observed over 40 clusters with the centimeter-wave 
SZ system during the 1994-1998 observing seasons. Only 
some of these clusters were observed for a significant 
amount of time; some observations were intended to survey 
for point sources and to define a sample for future work. 
To date, over 25 clusters have been detected significantly; 
analysis of a sample of 18 are presented here. 

The cluster targets were selected from a flux-limited, ho- 
mogeneous sample of clusters (Stocke et al. 1991; Gioia & 
Luppino 1994; Nichol et al. 1997) identified from the Ein- 
stein Extended Medium Sensitivity Survey (EMSS) (Gioia 
et al. 1990) and from two flux-limited samples (XBACS, 
Ebeling et al. (1996); BGS, Ebeling et al. (1998)) from 
the ROSAT All-Sky Survey, as well as public ROSAT 
data. Identifying clusters based on X-ray emission rather 
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than galaxy surface-density enhancements ameliorates the 
problems of false detections due to chance superpositions 
and of missed clusters due to smaller-than-average back- 
grounds. 

We selected massive clusters for our sample. X-ray stud- 
ies of clusters in David et al. (1995) and Mohr et al. (1999) 
indicate that the gas mass fraction near the virial radius 
increases as cluster mass increases, but that above ~ 5 
keV the gas mass fraction at the virial radius is constant. 
At the initiation of this work, X-ray temperatures were 
not widely available for distant clusters. We chose instead 
to select clusters on the basis of luminosity. We expect 
a cluster's X-ray luminosity to be a better predictor of 
mass than X-ray surface brightness, as it will be less sen- 
sitive to projection effects and contamination by cooling 
flows and dynamical activity in the ICM. Although cooling 
flows have been observed to contribute as much as 70% of a 
cluster's luminosity, typically they only contribute 10-30% 
(Peres et al. 1998). Subsequent X-ray spectral measure- 
ments conflrm that the clusters in this sample all have 
emission-weighted temperatures greater than 5 keV and 
therefore qualify as massive clusters for our purposes. 

Our SZ observing scheme requires the clusters to be at 
declination greater than —10°. The apparent size of the 
cluster must also be small enough so that the angular size 
is comparable to the spatial frequencies the interferometer 
samples; this is generally satisfied if the cluster redshift 
is greater than about 0.14. For the initial cluster obser- 
vations, we did not pursue observations towards cluster 
fields which hosted point sources with flux densities greater 
than ~10 mJy; fewer than ~15% of cluster fields had such 
point sources. We have since confirmed that we can reli- 
ably remove such point sources from the data, and we are 
pursuing observations towards these fields. 

It is possible that selecting against clusters with strong 
point sources may introduce a bias. This bias would be 
redshift dependent because, while the SZ effect magnitude 
is not diminished by distance, the flux of a point source 
associated with the cluster is. Clusters with radio-loud 
central galaxies will be less likely to be dropped because 
of point source contamination if they are distant. Peres 
et al. (1998) study a sample of 55 nearby X-ray clusters, 
40% with inferred cooling flow mass deposit rates of over 
100 Mq/yi. Forty-one of these clusters have radio detec- 
tions or upper limits at 1.4 or 5 GHz, and 33 of these have 
detected radio flux. Peres et al. (1998) cross-correlate the 
radio data and find only a weak correlation between the ra- 
dio power of the brightest cluster galaxy and the strength 
of the cooling flow. They do flnd that the largest cooling 
flows have the strongest radio fluxes, though. By selecting 
against clusters with very strong radio emission, we may 
be removing from our sample clusters which have not un- 



dergone recent mergers strong enough to disturb a cooling 
flow. Again, we expect this effect to be small, if present, 
as 85% of clusters were kept in the sample. 

The 18 clusters in our sample are listed in Table 1, along 
with the published redshifts and we used in the fg anal- 
ysis, and the X-ray luminosities. SZ images of the clusters 
are presented in Figure 1. ordered by redshift. We note 
that the quality of the detection reflects the r.m.s. sensi- 
tivity of the observation and the intrinsic strength of the 
SZ effect and not the cluster's distance. A Gaussian taper 
is applied to the u — v data to emphasize the structure 
on large scales. This taper depends on the range of u-v 
radii in each cluster's data set; the tapers are generally 
0.9 to 1.2 kA. Higher resolution images can be made from 
these data in order to emphasize smaller cluster struc- 
tures, e.i?., Carlstrom et al. (1996). Because the primary 
beam for the BIMA system is considerably larger than that 
for the OVRO system (396" and 252", respectively), the 
images produced from BIMA data include information on 
the decrement at larger scales. The images are plotted in 
contours of 1.5 a, and the restoring beam is shown in the 
bottom left corner of each image. 

The interferometric SZ data is necessarily spatially fll- 
tered; the visibility function will not be measured at every 
spatial frequency. The images are used to indicate the 
signal-to-noise ratio of the cluster detections, but we do 
not fit models to the images. 

3. SZ GAS MASS AND GAS MASS FRACTION 
MEASUREMENT METHODS 

3.1. Model 

We compare a model to the data in the spatial frequency 
domain, where the noise characteristics and the spatial fil- 
tering of the interferometer are well-understood. 

We fit to a /?-model (Cavaliere & Fusco-Femiano 1976, 
1978), which has been widely used to fit the density and 
temperature profiles of cluster galaxies and the ICM. We 
make the simplifying assumptions that the cluster gas is 
isothermal and the density distribution is spherically sym- 
metric. We consider the effects of these simplifications on 
our results in Section 3.5. 

In this model, the electron number density as a function 
of radius, r, takes the form: 

/ ^2 X -3/3/2 
rieir) = Heo U + ^ j , (4) 

where rc, the core radius, and /3 are fit parameters, and 
Ueo is the central electron number density. For isothermal 
gas with temperature Tg, Equation 4 predicts the following 
two-dimensional SZ temperature decrement: 

1 3/3 

AT{9)^AT{Q)i^l + -j , (5) 

where 9 = t/Da, Da is the angular diameter distance, 
6c = Tel Da, and AT(0) is the temperature decrement at 
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zero projected radius. The central electron density can 
therefore be recovered from this relation: 



-AT(0)mec2 1 
TcMB 2kaT Te 



+00 / 



'30/2 



(6) 

where the integral, dl, is along the line of sight. The mean 
molecular weight is assumed to be constant througout the 
gas, so the electron number density, rig, should trace the 
gas density. 

3.2. Fitting Procedure 

The change in spectral intensity of the CMB due to the 
Sunyaev-Zel'dovich effect is calculated for the Rayleigh- 
Jeans approximation (c./. Rephaeli (1995); Challinor & 

Lasenby (1998)): 

ATsz 



TcMB 



2 T 

yx e 



[xcoth{x/2)-A + eef{x)], (7) 



Rj 

where x = hu/kTcMB and 6e = kTf,/me(?- We adopt the 
COBE FIRAS value of Tcmb = 2.728 K (Fixsen et al. 
1996). The last term, 9ef{x), corrects for relativistic ef- 
fects. At 28.5 GHz, ATsz/TcMB = -1-92 y in the non- 
relativistic Rayleigh-Jeans approximation. Including the 
relativistic correction for a temperature typical of massive 
clusters, kTe = 7 keV, ATsz/Tcmb = -1-88 y. 

The data are components of the Fourier transform of the 
sky brightness distribution, i.e., a measured amplitude and 
phase for each two-dimensional spatial frequency, or u-v 
pair, sampled. The model is constructed in image space 
by filling out a regular grid with the SZ model (Equa- 
tion 5) multiplied by the primary beam. This SZ image is 
fast Fourier transformed and the model is interpolated to 
the u-v data points to compare with the data using the 
statistic. The cluster center, /3, 6c, and Ar(0) arc allowed 
to float to find the minimum using a downhill simplex 
(Press et al. 1992). The position and flux density of any 
radio-bright point sources are also fitted. Since the pri- 
mary beam attenuation at any given point differs between 
the OVRO and BIMA datasets, and the intrinsic point 
source flux can vary with time, the point source fluxes for 
each dataset are allowed to vary individually. 

The fits are performed jointly on all datasets for a given 
cluster. The shortest telescope spacing corresponds to the 
shadowing limit; for OVRO data this limit is 1 kX, for 
BIMA data this is 0.58 A;A. We use the holographically de- 
termined primary beam models when modeling the data, 
and the entire datasets are used to do the analysis. 

3.3. Constraints on Fit Parameters 

The cluster's centroid position and the point source 
fiuxes and positions are well constrained by the data. The 
fitting program consistently obtains the same values for 



the centroid positions. The initial guesses for the point 
source parameters are made using DIFMAP (Pearson et al. 
(1994)), an interactive mapping program, to inspect the 
high spatial frequency (|u^ + u^j > 2.0 kA) data in which 
the SZ effect contributes very little signal. The uncertainty 
introduced by point sources into the ICM parameters is 
discussed in Section 3.5.2. 

The cluster centroid and point source fiuxes and po- 
sitions are fixed to their best-fit values while the cluster 
shape parameters are fitted. We found no appreciable vari- 
ation of best fit centroid position with shape parameters. 
To illustrate the constraints these data place on (3 and 
Oc a grid search is performed over these parameters with 
AT(0) allowed to assume its best fit value at each grid 
point. In Figure 2, we show the confidence intervals for 
P and 9c for a representative cluster, Abell 1995. The 
solid contours indicate Ax^ = 2.3, 4.61, and 6.17 which 
enclose regions corresponding to 68.3%, 90.0%, and 95.4% 
confidence, respectively, for the two-parameter fit. The 
projection of the dashed lines, Ax^ = 1.0, 2.71, and 6.63, 
indicate the 68.3%, 90% and 99% confidence interval on 
the single parameter. At each (/?, 6c) point, the width of 
the 68% confidence interval for Ar(0) is about 10-15% of 
the best fit Ar(0) value. In Patel et al. (2000), we fit the 
ROSAT HRI data for this cluster, and find the fit values 
to be consistent with the SZ effect values. 

From the figure it is clear that /? and 6c are strongly 
correlated and the fit parameters (3 and 6c are not well- 
constrained individually by the SZ eff'ect data. 

3.4. Gas Mass Fraction Measurements 

The number of electrons in a given volume can be calcu- 
lated by integrating Equation 4. To recover the ICM mass, 
we multiply by the proton mass and the nucleon/electron 
ratio of 1.16. To extract the central electron density, Uco, 
for a given set of model parameters and measured electron 
temperature, we perform the integral in Equation 6. For- 
mally, this integral extends from the observer along the 
line of sight through the cluster infinitely; in practice, a 
cutoff radius of 8-10 cluster core radii is used. 

We note that although the fit parameters /3, 6c and 
AT(0) are not constrained strongly individually, the com- 
bination of these three parameters does constrain the gas 
mass quite well. This follows from the fact that the SZ ef- 
fect is, under the isothermal assumption, a direct measure 
of the gas mass on the scales to which our observations arc 
sensitive. We present gas masses for the 18 clusters in our 
sample in Table 2. To convert angular sizes to lengths, we 
have assumed h = 0.7, Om= 0.3, Oa=0.7. 

The distribution of the cluster's total mass, mainly com- 
prised of dark matter, can be inferred from the modeled 
gas pressure distribution, since the temperature of the gas 
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and its spatial distribution are constrained by the cluster's 
gravitational potential. We make the assumption that the 
gas is in hydrostatic equilibrium in this potential and that 
bulk flows and other non-thermal processes do not con- 
tribute significantly to the gas pressure. Under the as- 
sumptions of spherical symmetry and isothermal gas, the 
total mass of a cluster within radius r, is 



M(r) 



3kTeP 



(8) 



Gfinip + 

where /imp is the mean molecular weight of the gas. To 
calculate fi, we assume the gas has solar metallicity as 
measured by Anders & Grevesse (1989) and that fi is con- 
stant throughout the gas. The value of n will change 3-4 
% depending on the solar metallicity measurements one 
adopts; the metallicity in clusters is not well known, how- 
ever, and although an incorrect choice for /x will introduce 
a systematic error, it will be much smaller than the sta- 
tistical errors involved. Note from Equation 8 that the 
total mass depends only on the shape of the gas distri- 
bution, and is independent of the value of the central gas 
density, and therefore of the uncertainties in AT(0). Us- 
ing the derived shape parameters, (3, 9c, and the measured 
gas temperature, we derive the total mass, denoting it the 
"HSEmass". 

To measure the quantities of interest and their associ- 
ated uncertainties, we determine an appropriate range /3, 
6c, and AT(0) for each cluster with a coarse grid, and 
then construct a finer grid near the best fit parameter val- 
ues. The cluster's gas mass, HSE mass, gas mass fraction, 
and statistic are derived at each grid point. The 68% 
confidence interval on each quantity is determined from 
the range contained in the x^(best fit) - = ^X^< 1-0 
volume of the parameter grid. We prefer to measure the 
masses and mass fractions in the largest volume permitted 
by our method, since the fair sample assumption is best 
at large radii. The largest spatial scale on which we can 
constrain the model depends on the u-v points at which we 
significantly detect signal. To determine this scale, we cal- 
culate the statistical uncertainties in the fg measurement 
due to the shape parameter uncertainties for a number of 
radii from 10" to 150". We find we best constrain fg when 
it is calculated within a radius of around 65" (see Grego 
(1999)). The gas masses and gas mass fractions within 
a 65" radius along with their associated 68% confidence 
intervals are presented in Table 2. The gas mass and fg 
results depend on the assumed cosmology through the an- 
gular diameter distance, Da- For the gas mass fractions 
reported, we use J7m = 0.3, Qa = 0.7. 

The SZ gas mass is inversely proportional to the as- 
sumed electron temperature: Mgas{SZ) oc l/Tg and the 
HSE total mass measurement is directly proportional to 
Te- Mhse oc Tg. The gas mass fraction then is quite sen- 
sitive to temperature: fg oc l/T^ . The uncertainties from 



the temperature measurement are of the same order as the 
statistical uncertainties from the SZ model fitting at the 
lower redshifts, and dominate the SZ uncertainties for the 
most distant clusters. 

3.5. Systematic Effects 

3.5.1. Emission- Weighted Temperature 

When available, we have used emission-weighted tem- 
peratures which were examined and corrected for the pres- 
ence of cooling flows. The central surface brightness excess 
exhibited by many clusters is interpreted as emission from 
centrally-concentrated dense gas, e.g., Fabian (1994), the 
cooling time of which is shorter than the Hubble time. 
Such cooling flows can bias the emission-weighted temper- 
atures lower than the density-weighted or virial tempera- 
ture of the cluster. Allen & Fabian (1998) find that mod- 
eling clusters with a cooling flow spectral component in 
addition to the thermal component significantly reduces 
the scatter in the luminosity-temperature relation. We 
have used these cooling flow-corrected temperatures where 
available. 

The emission-weighted ICM temperatures we have 
adopted from the literature may also have errors due to 
contamination from other sources in the field. The ASCA 
observatory was the source for most of the published ICM 
temperatures we use in this work. As its half power diam- 
eter is ~ 3', it is nearly impossible to remove the effects 
of point sources on spectra of distant clusters obtained 
with ASCA. Since the measurement is so strongly depen- 
dent on an accurate measurement of Tg, this is likely to be 
the largest source of systematic uncertainty. Fortunately, 
many of the clusters in our sample are scheduled to be ob- 
served with the Chandra and XMM observatories, which 
will be better able to distinguish ICM emission from point 
source emission and toconstrain the ICM temperatures. 

3.5.2. Radio Point Sources 

We detect radio-bright point sources in about half of the 
observed clusters. The point sources with fluxes exceeding 
three times the r.m.s. of the high resolution (> 2000 A) 
maps can be reliably identifled from the SZ data. We esti- 
mate the maximum effect of undetected point sources by 
adding an on-center point source to a representative clus- 
ter data set and fitting this new data set not accounting 
for the added point source. We place a 3 cr point source 
at the cluster center where typical r.m.s. sensitivities in 
the OVRO and BIMA high resolution maps are roughly 
61 /xJy and 163 /xJy respectively. This point source causes 
the magnitude of the decrement to be underestimated (and 
thus the gas mass fraction too) by 15% for OVRO data and 
20% for BIMA data. Such a point source at the cluster 
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center is highly unhkely but places limits on the maximum 
effect from undetected point sources. 

3.5.3. Departures from an Isothermal, Spherical ICM 

Our assumptions that the intra-cluster medium is 
isothermal and spherical are at some level approximations. 
Markevitch et al. (1998) report moderate temperature gra- 
dients in a sample of 30 nearby clusters, although in a sim- 
ilar analysis, Irwin ct al. (1999) do not find such structure. 
Neglecting to account for existing temperature gradients 
in the ICM may systematically affect the gas and HSE 
masses. If such a gradient is present, the true tempera- 
ture in the central region may be higher than the emission- 
weighted temperature we use, and the fitted shape param- 
eters from the isothermal SZ analysis may no longer accu- 
rately describe the density distribution. As yet, there are 
no strong observational constraints on temperature struc- 
ture in clusters beyond z = 0.1, as there have been no suit- 
able telescope facilities for the task. We anticipate that 
the Chandra and XMM X-ray observatories will greatly 
improve this situation. 

Our observation scheme provides information on the 
two-dimensional decrement, and we observe that the clus- 
ters are not strictly spherical. For this sample of clusters, 
we find the mean of the best-fit axis ratios to be 0.89 ± 
0.12. In previous work (Grego et al. (2000); Grego (1999)), 
we relaxed this assumption and permitted the density dis- 
tribution to be cUipsoidally symmetric, but the unknown 
orientation and three-dimensional geometry introduce a 
large uncertainty in the HSE mass. For a sufficiently large 
sample chosen without orientation bias, deviations from 
spherical symmetry will not strongly affect the results. As 
a point of comparison, the effects of cluster shapes on de- 
terminations of the cluster size have been investigated in 
Sulkanen (1999). He calculates Hubblc's constant in a 
sample of simulated tri-axial clusters by comparing their 
predicted SZ and X-ray images. The X-ray flux from the 
cluster at any point in the sky is proportional to / nidi, 
integrated along the line of sight through the cluster, while 
the SZ effect is proportional to / n^dl, so the two can be 
compared to derive the size scale of the cluster; when this 
size scale is compared to the apparent size on the sky, the 
cluster's distance and hence Ho can be inferred. Sulkanen 
finds that when the images are fit by an spherical beta 
model with a core radius equal to the arithmetic mean of 
the two core radii from an elliptical fit, the recovered Ho 
for a sample of 25 clusters is unbiased. 

In an ongoing analysis of an ensemble of hydrodynami- 
cal cluster simulations, we also find that we do not intro- 
duce serious error with these assumptions. These sim- 
ulated clusters are produced within both low and high 
f^M cosmological models, and the temperature and den- 



sity structure is appropriate for cluster populations expe- 
riencing merging similar to that observed at redshifts z< 
0.1 (Mohr et al. 1995). We produce mock BIMA obser- 
vations of simulated clusters at the redshifts z=0.2 and 
z=0.6. Isothermal, spherical /3-model analyses of these 
SZE observations produce unbiased estimates of the ICM 
mass enclosed within the radius r2ooo, which roughly cor- 
responds to the scales measured in this experiment (Mohr 
et al. 2000). Should temperature and density structure in 
distant clusters be similar to that in the local sample, it 
should not be a significant source of systematic uncertainty 
or bias in our measurements. 

3.5.4. Validity of the HSE Approximation and 
Non-thermal Pressure Support 

Our method of measuring the total mass assumes the 
ICM is in hydrostatic equilibrium in the cluster poten- 
tial and supported only by thermal pressure. One test 
of this assumption is to compare the HSE-derived total 
cluster mass to the total mass derived from gravitational 
lensing models. Some mass comparisons in the literature 
(Miralda-Escude & Babul (1995); Loeb & Mao (1994); Wu 
& Fang (1997)) have suggested that the HSE method may 
systematically underpredict the cluster's total mass by a 
factor of 2-3, compared to a strong gravitational lensing 
analysis. Suggested explanations for this discrepancy in- 
clude elongation of the cluster along the line of sight and 
temperature structure in the ICM, which we discuss in 
Sections 3.5.3, and non-thermal pressure support of the 
gas in the cluster core. 

Further work has suggested that the details of the anal- 
ysis can have a significant effect, and may resolve the 
discrepancy. A weak lensing analysis was performed by 
Squires et al. (1997) on the cluster Abell 2218, which ap- 
pears to have discrepant masses in each of the three anal- 
yses above. This analysis, at larger radius than the strong 
lensing analyses, show the two methods predict masses 
which are consistent within the experimental uncertain- 
ties. In an examination of a sample of 13 clusters, Allen 
& Fabian (1998) finds that the lensing and HSE masses 
agree for clusters which appear to have a strong central 
cooling flow, when the cooling flow is taken into account; 
in these clusters, the X-ray and lensing core radii are con- 
sistent with each other, and the mass agreements suggest 
HSE is a reasonable approximation. For clusters with- 
out strong cooling flows, the X-ray core radii arc gener- 
ally larger than the lensing radii, and offsets between the 
centers of the distributions are observed, suggesting that 
HSE is not appropriate for the cluster cores. Outside the 
cluster core (at radii 400 kpc, weak lensing and X- 
ray masses are consistent with each other both for cooling 
flow and non-cooling flow clusters. Lewis et al. (1999) 
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compare X-ray HSE masses to the dynamical masses cal- 
culated from the galaxy velocity dispersions, and find the 
average Mdyn/^HSE to be 1.04±0.07, which also suggests 
the HSE method does not introduce a systematic bias. 

Possible sources of non-thermal pressure in the ICM are 
bulk flows and magnetic fields in the gas. Intracluster 
magnetic fields are typically a few /iG (Kim ct al. (1990); 
Taylor & Perley (1993)), an order of magnitude smaller 
than the level at which the fields would contribute signif- 
icantly to the dynamics of cluster gas, although stronger 
fields, ^ 10—100 /iG, have been measured in a few clus- 
ters (Taylor & Perley (1993)). There is some evidence for 
the persistence of bulk flows in clusters undergoing merger 
events (Bliton ct al. (1998)). It remains to further in- 
vestigation how significant a role these effects play in the 
physics of cluster gas, but currently there is no evidence to 
suggest a significant systematic error in the HSE method. 

3.5.5. Cosmic Microwave Background Anisotropics 

The SZE observations arc also sensitive to intrinsic 
and secondary anisotropy in the Cosmic Microwave Back- 
ground (CMB) radiation. The theoretical expected and 
observed level of CMB anisotropy at the angular scales 
corresponding to those used for the SZE measurements 
presented here is small and safely ignored. The contribu- 
tion of primary anisotropy for a window function appropri- 
ate for our shortest baselines, for which the contribution 
would be strongest, has been calculated by Holzapfel et al. 
(2000a). For a flat universe, as indicated by recent CMB 
observations (Miller et al. 1999; de Bernardis et al. 2000; 
Hanany et al. 2000) the r.m.s. temperature fluctuations 
contributed by primary CMB anisotropy within our maps 
should be of order 2/iK or less. 

At the angular scales of our SZE measurements, sec- 
ondary CMB anisotropy is expected to be stronger than in- 
trinsic anisotropy. Holzapfel et al. (2000a) have tabulated 
the expected range of the magnitude of the temperature 
anisotropy due to the Visniac Effect and inhomogenous 
rcionization. The upper range for the r.m.s. temperature 
fluctuations at angular scales appropriate for the SZE mea- 
surements is only 5.6 /zK and 3.9 /xK, respectively. Added 
in quadrature, this gives an expected upper limit to the 
r.m.s. temperature fluctuations in our SZE observations 
of only 6.8 jjK. If present, this signal would contribute to 
the SZE maps as noise {i.e., would not lead to a bias in 
our SZE measurements). This level is muc;h smaller than 
the noise level obtained in our SZE observations. 

The dominant contribution to secondary anisotropy at 
the relevant angular scales is likely to be the SZE from 
undetected low mass clusters. Holder & Carlstrom (1999) 
estimate r.m.s. temperature fluctuations of order 2 /zK to 
12 /xK for the range of models they consider. Again, this 



level is small compared to our uncertainties, although ap- 
proaching the noise level in our deepest fields. 

It is unnecessary to depend on theoretical estimates of 
contributions from CMB anisotropy as we have direct mea- 
surements of 'blank' fields obtained with the same instru- 
ment as for the SZE observations (Holzapfel et al. 2000a). 
The r.m.s. level obtained on the deepest fields ranges from 
16 jiK to 20 /xK, just slightly above that expected from the 
instrumental noise. The most likely level of anisotropy, in- 
cluding undetected point sources, derived from the blank 
field data is 12 /uK and the 95% confidence level upper 
limit is 19 /xK. 

We conclude that temperature fluctuations due to pri- 
mary and secondary CMB anisotropy should have a negli- 
gible effect on the results derived from the SZE measure- 
ments reported here. 

4. RESULTS AND DISCUSSION 

4.1. Gas Mass Fractions 

As we discussed in Section 3.4, we measure the gas mass 
fraction within a fixed angular radius, which results in the 
measurements being made at different physical scales for 
clusters at different distances. To compare the gas mass 
fractions of different clusters, and to derive a result useful 
for cosmological tests, we scale the results for each clus- 
ter to a fiducial radius. An analytical expression for the 
variation in fg with radius is suggested by Evrard (1997), 
based on results in Evrard et al. (1996) and found to be 
consistent with the fg variation reported in the David et al. 
(1995) sample. We use a modified version of this to ex- 
trapolate the gas mass fractions we measure at 65" to the 
gas mass fraction expected at rsoo, the radius at which the 
cluster's total mass density is 500 times the critical mass 
density, where the cluster's baryon fraction should closely 
refiect the universal value. The physical radius at which 
the overdensity is 500 depends on the cluster's temper- 
ature (a mass indicator), and also its redshift, since the 
critical density will change with z. The scaling expression 
is as follows: 

f,{rMTe))=Mrx)(^^^^y, (9) 

where ri = 0.17, fg{r5oo{Te)) is the gas mass fraction within 
r5oo, and rx is the radius within which the gas mass frac- 
tion is measured. We modify Evrard's expression for rsoo, 
derived for low redshift clusters, to include the change in 
the value of Pc with redshift; pc{z) = Pc{z = 0){H/Ho)^, 
where = [(1 + z)^n m + {1+z)^{1-Qm -i^A) + ^A] 
{c.f. Peebles 1993, Eqn. 13.3): 

rM = (1.24 ± 0.09) [j^^^^^, ) h-^M,o. 

(10) 

The gas mass fractions within rsoo as derived by this 
relation are presented in Table 2. Figures 3a. and b. 
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show the gas mass fractions at r^oo as a function of and 
redshift. We see no correlation of gas mass fraction with 
temperature. We see no significant evolution of fg with 
redshift. Since fg depends on the cluster distance, fg oc 
T>A, and therefore the chosen cosmology, measurement of 
the gas mass fraction over a range of redshifts could be 
used in principle to constrain cosmological models. 

We calculate the mean gas mass fraction for the en- 
tire cluster sample, and derive the 68% confidence interval 
from the Ax^ statistic to a constant-value fit. For the en- 
tire sample, assuming r2M=0.3, r^A=0.7, we find the mean 
gas mass fraction to be 0.08lt^;^??/i^oo- 

We also calculate 

the mean and uncertainty for fg in the full sample, us- 
ing two alternative cosmologies, (riM=0.3, f2A=0.0) and 
(17m = 1-0, riA=0.0), to calculate the distances and scaling 
relation. In Table 3, we report these and the associated re- 
duced chi-squared (Xrerf) statistics, which range from 1.021 
to 1.056 for the full sample fits. The xled values do not dif- 
fer significantly enough to discriminate between cosmolo- 
gies, and it is clear that currently the uncertainties are too 
large for a cosmological test via geometry. 

We also calculate the mean fg in a homogeneous sub- 
sample of five clusters. These clusters are the five most 
luminous clusters in the flux-limited EMSS cluster sam- 
ple with z > 0.26 and declination > -10°: MS0451.6- 
0305, MS1137.5-H6625, CL0016-hl6, MS1358.4-h6245, & 
MS1054. 4-0321. For our standard cosmology, we find the 
mean in this sample to be 0.089lg;Q^g/ij~Qo. In all three cos- 
mologies, the gas mass fraction in the homogenous sample 
is consistent with the full sample value. 

We compare these SZ-derived gas mass fractions to other 
SZ-derived fg measurements. Recent cluster gas mass frac- 
tion measurements from SZ effect observations are pre- 
sented in Myers et al. (1997). In this work, the integrated 
SZ effect is measured using a single radio dish operating 
at centimeter wavelengths. The integrated SZ effect is 
used to normalize a model for the gas density from pub- 
lished X-ray analyses, and this gas mass is compared to 
the published total masses to determine the gas mass frac- 
tion. For three nearby clusters, A2142, A2256 and the 
Coma cluster, Myers et al. (1997) find a gas mass frac- 
tion of (0.061 ± 0.011)/iioo at radii of 1-1.5 h^Q^ Mpc; for 
the cluster Abell 478, they report a gas mass fraction of 
(0.16 ± 0.014)/ifo^o. 

4.2. Comparison of SZ and X-ray Results 

Gas mass fractions derived from X-ray images for a 
large, homogeneous, nearby sample of clusters are pre- 
sented in Mohr et al. (1999). For a subsample of 28 clus- 
ters with Te > 5 keV, they find the mean gas mass fraction 
within r5oo to be (0.0749±0.0021)/i^oo/^ at 90% confidence. 
The gas mass fractions derived from SZ measurements de- 



pend differently on the cosmology assumed than those de- 
rived from X-ray images, and this should be noted when 
comparing the results. 

Qualitatively, though, the comparison docs not suggest 
any large systematic offsets. This is a significant result, 
because a large clumping factor, c =< >^/^ / < rig >, 
has been suggested as an explanation for the high gas mass 
fractions in clusters (White et al. 1993; Evrard 1997). A 
cluster with clumping factor c would only require 1/c as 
much gas mass to produce the observed emission, and so 
the SZ and X-ray gas mass fraction measurements would 
differ by a factor of ~ c. 



4.3. Comparison of Baryon Fraction with CIb 

The relative abundance of deuterium and hydrogen pro- 
vides a particularly strong constraint on the baryonic mat- 
ter density (Copi et al. 1995). A firm upper limit to 
is set by the presence of deuterium in the local inter- 
stellar medium. This constrains the value of Qb to be 
less than 0.031/ij^oo (Linsky et al. 1995). Measurements 
of the D/H ratio in metal-poor Lyman-a absorption line 
systems in high-redshift quasars put an even more strin- 
gent constraint on the baryonic mass density. For this 
analysis, we adopt the published value at 95% confidence 
from the Buries & Tytler (1998) absorption line anaylsis, 
Os = (0.019 ± 0.002)/ifcfo. 

We can use the gas mass fractions to find the value of 
Qm in a self-consistent manner. In Figure 4, we show the 
value of fijvf implied by the measured gas mass fractions 
when we assume a flat universe {flA= 1- J^m) and h = 0.7 
to calcinate the angular diameter distance and r^oo scaling 
factor from Equation 9: < ^b/ fs/hro . The upper 
limit to Om and its associated 68% confidence interval is 
shown as a function of Qm- The measured gas mass frac- 
tions are consistent with a flat universe and h = 0.7 when 
Om is less than 0.40, at 68% confldence. For our measure- 
ments to be consistent with Om = 1.0 in a flat universe, 
the Hubble constant must be very low, h less than ^ 0.30. 

For a more realistic estimate, we could include the 
baryon contribution from galaxies, and attempt to ac- 
count for the overall dimunition of the baryon fraction in 
clusters with respect to the universal value, since some 
baryons are expected to not become bound to the clus- 
ter. Following White et al. (1993), we estimate the galaxy 
mass to be a fixed fraction of the cluster gas, with the 
same fraction as is observed in the Coma cluster, Mb = 
Mgasi^ + 0.20/i^/^). For a realistic equation of state, the 
gas in the cluster will be more extended than the dark 
matter and the baryon fraction at rsoo will be a modest 
underestimate of the true baryon fraction (Evrard 1997), 
fg {^500) = 0.85 X /(.(universal). These assumptions lead 
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to 

/s = (/,(l + 0.2/i3/2)/o.85). (11) 

Using this estimate of the baryon fraction, and h = 0.7 in 

a flat cosmology, in Figure 4 we show our best estimate 
of ilAf as a function of cosmology. Thus we find our best 
estimate of is ~0.25. 

4.4. Future Work 

There are several improvements to this work which will 
be made in the near future. More clusters will be added to 
the sample as SZ observations continue. And the potential 
also exists for improving the centimeter-wave SZ interfer- 
ometer system dramatically by taking advantage of the 10 
GHz output of the SZ receivers; currently a maximum of 
2 GHz are correlated at OVRO and effective bandwidth of 
0.5 GHz are correlated at BIMA. 

One of the main sources of uncertainty in these measure- 
ments originates in the emission- weighted gas temperature 
measurements; as fg oc T~^, the 10-20% uncertainties in 
Te roughly double to 20-40% uncertainties in the gas mass 
fraction. A large number of these clusters are scheduled 
to be observed in the Chandra X-Ray Observatory GTO 
and GO phases, which should improve the situation con- 
siderably. 

Numerical simulations will also help identify other 
sources of systematic error incurred in the observational 
and analysis program. An analysis is in preparation of hy- 
drodynamic simulations of a sample of clusters to quantify 
any biases we may introduce to the gas mass fraction mea- 
surements with the interferometric method and through 
the assumptions we make in the fitting and analysis of the 
clusters. 
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Table 1 
The Cluster Sample 
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References. — A92 Allen (1992); AF Allen & Fabian (1998); AE Arnaud & Evrard (1998); 
AEest estimated from L^-T relation of AE; BA Bade, N. et al. (1998); C95 Crawford et al. 
(1995); D98 Donahue et al. (1998); D99 Donahue et al. (1999): EB Ebcling et al. (1996), E98 
Ebeling (1998); GI Gioia et al. (1990); HB Hughes & Birkinshaw (1998); LB LeBorgne et al. 
(1992); LG Luppino & Gioia (1995); M88 Mellier et al. (1988); MS Mushotzky & Scharf (1997); 
OT Ota et al. (1998); PA Patel et al. (2000); SR Struble & Roo d (1991) 
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Table 2 

SZ-DERivED Gas Masses and Mass Fractions, using = 0.3, = 0.7 
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Table 3 
Mean Gas Mass Fractions 



nM=o.'d,nA=o.7 nM=o.3,nA=o.o riM=i-o,fiA=o.o 



sample number of clusters /g/iioo Xred fghioo xled fghwo xled 



full sample 18 0.0811^;^^^^ 1.021 0.074^^;^^^ 1.027 0.0681^;^^^ 1.056 

EMSS subsample 6 0.089^^19 1-208 0.077^^16 1-258 0.067^^14 1-352 



Note. — The mean gas mass fractions for the noted samples are presented. The gas mass fractions 
depend on DJ^ and so the results are presented for three sets of cosmological parameters. {Ho = lOO/i 
km s~^ Mpc~^). 
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10 20 30 40 50 60 70 80 90 100 110 120 

0^ (arc seconds) 

Fig. 2. — Confidence intervals for /3 and 9c from a /3-model fit to tiie Abell 1995 data. At eacli point in tfie plot, the central decrement was 
allowed to assume its best fit value. The solid contours arc marked for Ax^ = 2.3, 4.61, and 6.17 which indicate 68.3%, 90.0%, and 95.4% 
confidence, respectively, for the two-parameter fit. The dashed lines indicate Ax'^ =1.0, 2.71, and 6.63. The projection onto the /? or 6c ctxis 
of the interval contained by these contours indicate the 68.3%, 90% and 99% confidence interval on the single parameter. 
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Fig. 3. — Gas mass fractions within rsoo for the entire sample, assuming the cosmology {Qm=0.3, f2A=0.7). a) Gas mass fraction vs. ICM 
emission-weighted temperature, b) Gas mass fraction versus redshift. 
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Fig. 4. — Upper limit on the total matter density, iiM< f2s/(/s/i7o) (full line) and its associated 68% confidence region (dotted lines), 
as a function of cosmology {Qa= 1 - ^m)- The baryon fraction is a function of cosmology through the angular diameter distance relation 
as well as the scaling relation in Equation 9. The intercept between the upper dotted line and the dashed line [(^m= f2s/(/B'i7o)] gives the 
upper limit to Qm at 68% confidence. The dot-dashed line shows the total matter density when the baryon fraction includes an estimate of 
the contribution from baryons in galaxies and those lost during cluster formation. The intercept of the dot-dashed line and the dashed line 
gives our best estimate of flu (~0.25) assuming a flat universe with h=0.7. 
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